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The production of single photons in Pb + Pb collisions 
at the CERN SPS as measured by the WA98 experiment is 
analysed. A very good description of the data is obtained if a 
quark gluon plasma is assumed to be formed initially, which 
expands, cools, hadronizes, and undergoes freeze-out. A rich 
hadronic equation of state is used and the transverse expan- 
sion of the interacting system is taken into account. The re- 
cent estimates of photon production in quark-matter (at two 
loop level) along with the dominant reactions in the hadronic 
matter leading to photons are used. About 50% of the single 
photons are seen to have a thermal origin. An addition of the 
thermal and prompt photons is seen to provide a very good 
description of the data. Most of the thermal photons having 
large transverse momenta arise from the quark-matter, which 
contributes dominantly through the mechanism of annihila- 
tion of quarks with scattering, and which in turn is possible 
only in a hot and dense plasma of quarks and gluons. The 
results thus confirm the formation of quark gluon plasma and 
the existence of this mechanism of the production of single 
photons. 

PACS numbers: 25.75.-q, 12.38Mh 



The search for quark-gluon plasma, which filled the 
early universe microseconds after the big bang and which 
may be present in the core of neutron stars, is one of 
the most notable collective efforts of the present day nu- 
clear physics community. Its discovery will provide an 
important confirmation of the predictions of the statisti- 
cal Quantum Chromodynamics (QCD) based on lattice 
calculations. It has been recognised for a long time |l]] 
that electromagnetic radiations from relativistic heavy 
ion collisions in these experiments would be a definitive 
signature of the formation of a hot and dense plasma of 
quarks and gluons, consequent to a quark-hadron phase 
transition 0j . Once other signs of the quark-hadron tran- 
sition, e.g., an enhanced production of strangeness, a sup- 
pression of J/tp production, radiation of dileptons, etc., 
started to emerge Q , it was imperative that the more di- 
rect, yet much more difficult to isolate, signature of the 
hot and dense quark-gluon plasma, the single photons 
were identified. The WA98 experiment || has now re- 
ported observation of single photons in central Pb + Pb 
collisions at the CERN SPS. 

In the present work we show that these data are very 
well described if we assume that a quark-gluon plasma 
was formed in the collision. 

In order to put our findings in a proper perspective, 



let us recall that the publication of the upper limit of the 
production of single photons in S+Au collisions at CERN 
SPS B by the WA80 experiment was preceded and fol- 
lowed by several papers J5|,^| exploring their connection 
to the quark-hadron phase transition. An early work, 
by the present authors ||, reported that the data were 
consistent with a scenario where a quark gluon plasma 
was formed at an initial time To ~ 1 fm/c, which ex- 
panded and cooled, got into a mixed phase of quarks, 
gluons, and hadrons, and ultimately underwent a freeze- 
out from a state of hadronic gas consisting of 7r, p, to, and 
•q mesons. On the other hand, when the initial state was 
assumed to consist of (the same) hadrons, the resulting 
large initial temperature led to a much larger production 
of single photons, in gross violation of the upper limit. 

A reanalysis of the WA80 data on single photons was 
reported recently j?J which incorporated two important 
developments in the field during the last few years, which 
are worth recalling. Firstly, it was realized that the 
hadronic equation of state must be generalized to include 
all of the hadrons || ( limited to M < 2.5 GeV, in prac- 
tice). This was prompted and supported by the success 
of the thermal models in describing particle production 
in these collisions. This implied that the hadrons were in 
chemical equilibrium || at least at the time of (chemical) 
freeze-out. These hydrodynamical calculations have been 
shown to provide a very good explanation of the pr spec- 
tra measured by the NA49 and NA44 experiments jl0| . 

Secondly, an evaluation of the rate of single photon 
production from the quark matter to the order of two- 
loops was reported recently by Aurenche et al pd| , ^2[ . 
This had two quite important results: (i) a substantial 
contribution of the bremsstrahlung (q q (g) — > qq (g) 7) 
process for all momenta in addition to the Compton 
{Q (<?) .9 ~~ * 9(9)7) P ms annihilation (qq — > 177) con- 
tributions included in the one-loop calculations available 
in the literature |f^,[l4||, and (ii) a large contribution by 
a new mechanism which corresponds to the annihilation 
of a quark (scattered from a quark or a gluon) by an 
anti-quark. These new rates were shown [ p"5[ to lead to 
a considerable enhancement of the production of single 
photons at SPS, RHIC, and LHC energies, if the initial 
state is approximated as an equilibrated plasma. 

It was also reported |7j] that when allowances were 
made for the above considerations, the WA80 upper limit 
was still consistent with a quark hadron phase transition, 
while a treatment without phase transition was untenable 
as it involved several hadrons/fm 3 , at the initial time. 



We add that there can be a production of high mo- 
mentum single photons during the pre-equilibrium phase, 
when treated within the parton cascade model |l6| ] , from 
the fragmentation of time-like quarks (q — > qj) pro- 
duced in (semi)hard multiple scatterings 

The rate for the production of hard photons evaluated 
to one loop order using the effective theory based on re- 
summation of hard thermal loops is given by jl3],[l4| : 
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where the constant c rs 0.23. The summation runs over 
the flavours of the quarks and ej is their electric charge in 
units of charge of the electron. The rate of production of 
photons due to the bremsstrahlung processes evaluated 
by Aurenche et al is given by: 
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and the exp ressions for Jt, Jl, and I(E, T) can be found 
in Ref. Q. 

And finally the dominant contribution of the qq an- 
nihilation with scattering obtained by Aurenche et al is 
given by: 
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Note that all the three contributions turn out to be es- 
sentially of the order aa s pj] |. It has been pointed out 
recently fl^J that the values of Jt and Jl given origi- 
nally by Aurenche et al |H| / are too large by a numerical 
factor of 4- We use the corrected values in the following. 

The estimate of prompt photons is taken from the work 
of Wong and Wang which employs the NLO pQCD 
along with the inclusion of the effects of intrinsic partonic 
momenta (< k T >= 0.9 GeV 2 ; see discussion later). 

We assume that a chemically and thermally equili- 
brated quark-gluon plasma is produced in such collisions 
at the time To (see later), and use the isentropy condi- 
tion || ; 
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to estimate the initial temperature, where At is the 
transverse area. 

We have taken the average particle rapidity density as 
750 for the 10% most central Pb + Pb collisions at the 
CERN SPS energy as measured in the experiment. We 



estimate the average number of participants for the cor- 
responding range of impact parameters (0 < b < 4.5 
fm) as about 380, compared to the maximum of 416 for 
a head-on collision. We thus use a mass number of 190 
to get the radius of the transverse area of the colliding 
system and neglect its deviations from azimuthal symme- 
try, for simplicity. As this deviation, measured in terms 
of the number of participants, is marginal (< 9%) we ex- 
pect the error involved to be small. We also recall that 
the azimuthal flow is minimal for central collisions. 

We take a = 42.257r 2 /90 for a plasma of massless 
quarks (u, d, and s) and gluons, where we have put the 
number of flavours as ~ 2.5 to account for the mass of the 
strange quarks. We now use Eq.([|) to estimate the (aver- 
age) initial temperature, with the additional assumption 
of a rapid thermalization pp| so that the formation time 
is decided by the uncertainty relation and tq = l/3Xo. 
This Tq is then used to get the (average) initial energy 
density. 

It is important to have a proper initial energy density 
profile as it affects the hydrodynamic developments by 
introducing additional gradients. We assume it to follow 
the so-called 'wounded-nucleon' distribution, which for 
central collision of identical nuclei leads to: 



e(r ,r) oc / p(y r 2 + z 2 ) dz 



(5) 



where p is the (Woods-Saxon) distribution of nucleons 
in a nucleus having a mass number of 190 and r is the 
transverse distance. This is prompted by the experimen- 
tal observation that transverse energy deposited in these 
collisions scales with the number of participants. The 
normalization in the above is determined from a numer- 
ical integration so that: 
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We further assume that the phase transition takes 
place at T — 180 MeV and the freeze-out takes place 
at 120 MeV. This value of the critical temperature is 
motivated by the recent lattice QCD results which give 
values of about 170 - 190 MeV Q, and the thermal 
model analyses of hadronic ratios which suggest that the 
chemical freeze-out in such collisions takes place at about 
170 MeV. ( A recent analysis by Becattini et al yields a 
value of 181.3 ± 10.3 MeV Q for the chemical freeze-out 
temperature.) The phase transition should necessarily 
take place at a higher temperature. 

The rates for the hadronic matter have been ob- 
tained fl3|| from a two loop approximation of the photon 
self energy using a model where tt — p interactions have 
been included. The contribution of the A\ resonance is 
also included according to the suggestions of Xiong et 
al p2j . The relevant hydrodynamic equations are solved 
using the procedure p3| discussed earlier and an integra- 
tion over history of evolution is performed Iq] . 



In Fig. 1 we show our results. The dashed curve gives 
the contribution of the quark-matter and the solid curve 
gives the sum of the contributions of the quark matter 
and the hadronic matter. The NLO pQCD estimates 
for prompt photons pp are also given. We see that the 
thermal photons contribute to about 50% of the total 
yield of the single photons and that the sum of thermal 
and prompt photons provides a very good description to 
the data. We also note that at higher transverse mo- 
menta most of the thermal photons have their origin in 
the quark matter. 

How sensitive are the results to the choice of our pa- 
rameters? In Fig. 2, we show our results where we vary 
the transition temperature by ± 20 MeV. It is seen that 
the results at higher kx (which have their origin in earlier 
times) remain unaltered, though the yield at the lowest 
transverse momenta increases with the decrease in Tq- 

The initial time tq affects the results much more 
strongly, as increasing it lowers the initial temperature 
(Eq.[|). In Fig. 3(a) we show our results for To = 0.20, 
0.40, 0.60, 0.80 and 1 fm/c, corresponding to T = 335, 
265, 232, 210, and 196 MeV. A comparison of this figure 
with Fig. 1 shows that the data clearly favour a large 
initial temperature (and early thermalization) . Recall 
that the hydrodynamic flow of all the systems (having 
same dN/dy ~ TqTq) are known to be nearly identical 
at later times and thus affect the hadronic data only 
marginally; see Fig. 3(b). 

A very important outcome of these results (Fig. 1) is 
that a very large part of the thermal component of the 
single photons is seen to have its origin in the quark- 
matter itself! Recall that the new (and dominant) mech- 
anism of the annihilation of quarks with scattering, sug- 
gested by Aurenche et al. , is operative only if a hot and 
dense plasma is formed (see the detailed discussion in the 
Appendix in Ref. (ll]]). Thus these results confirm the 
existence of this mechanism and the formation of quark 
gluon plasma in such collisions. 

Even though we realize that the creation of a hot (con- 
fined) hadronic matter in thermal and chemical equilib- 
rium within To ~ 0.20 fm/c, consequent to nuclear colli- 
sion is highly unlikely pjj , we estimate the initial tem- 
perature for such a system from Eq.@ for the hadronic 
equation of state used here, as more than 260 MeV, when 
the hadronic density would be w 10 hadrons/fm 3 pi. We 
consider this very unphysical and unlikely. A larger for- 
mation time will give a much lower initial temperature 
and fail to explain the data. 

How are we to understand the use of To = 1 /3T « 0.20 
fm/c here (see also |^j) against the canonical value of 1 
fm/c, employed often? Firstly, within the model used, 
this value is favoured by the data (Fig. 3). Secondly, if 
a larger value of t is used, then an allowance should be 
made to supplement the predictions with an appropriate 
pre-equilibrium contribution (see e.g., Ref. p5| ). Thirdly, 
we note that the matter at z = starts interacting by 
t = —R/j ~ —0.7 fm/c in the present case, when the 
two nuclei start touching. Thus by the lapse of r = 0.2 



fm/c, the matter there has been under interaction for a 
time ~ 1 fm/c, which may be enough for the formation 
of the plasma. 

Finally, a very important confirmation of our findings 
comes from the observations of Eskola et al p^j, that 
that a saturation of partons signaling a complete filling 
up of the transverse area by coloured quanta in collision 
of lead nuclei at SPS energies is indeed attained when the 
momentum transfer in partonic collisions is of the order 
of 1 GeV leading to a temperature ~ 300 MeV at To ~ 
0.2 fm/c. 

We add that the model developed here provides a very 
good description |10| of the intermediate mass dilepton 
excess measured by the NA50 group. 

Are we justified in making the assumption of a chem- 
ically equilibrated plasma, considering that indeed the 
predictions at the lower transverse momenta are close to 
the upper limits given by the experiment? This needs 
to be investigated (see Neumann et al ||) as also the 
effect of (likely) medium modification of hadron proper- 
ties. The neglect of the baryo-chemical potential for the 
QGP is perhaps justified as the net-baryon to hadron 
ratio is quite small (^7), especially in the region of the 
central rapidity. Finally, we may add that the photon 
rates used in these calculations are strictly valid only for 
a s <C 1 and that the consequences of considering higher 
loops remains to be seen. 

Before summarizing, let us return to the question of 
prompt photons. A detailed discussion on them is be- 
yond the scope of this paper, and the debate on the re- 
producibility of single photons data in fixed target pp(A) 
experiments is inconclusive. If we are to believe the re- 
sults of Wong and Wang |l8|], which we have employed, 
then the prompt photons contribute about half of the 
total yield in the present work, As mentioned earlier, 
these results are obtained by using NLO pQCD predic- 
tions along with the inclusion of intrinsic momenta of 
partons. 

The classic paper of Owens Q discusses the need to 
account for the intrinsic transverse momenta of partons. 
However, that work also talks of the need to introduce 
a cut-off in Q 2 , below which pQCD can not be applied 
and to avoid singularities in the parton-parton matrix 
elements. This discussion is absent in Ref. |l8|,|3(|. Large 
enhancements can be obtained depending on the the cut- 
offs employed and the < kr > used. The extent to which 
these considerations will affect the results of Wong and 
Wang is not known. 

We also recall the exhaustive work of Vogelsang and 
Whalley p3[ ; especially their Fig. 31, where the differ- 
ences of all the pp data from NLO pQCD predictions are 
plotted. It is seen that while all the pp data for y/s > 
23 GeV are quantitatively explained by NLO pQCD, the 
one at 19.4 GeV is underpredicted by a factor of 4-5. 
If this trend is assumed to continue then at the relevant 
nucleon-nucleon energy of 17.3 GeV, this difference would 
mount to a factor > 10 (!). The NLO pQCD analysis of 
these authors is also applied to p + Be data at 31 GeV 



(by normalizing it to pp) and a very good description is 
obtained, while the same procedure underestimates the 
p+C data at 19.4 GeV by a factor of 4-5. Several of these 
data have also been critically examined by Aurenche et 
al. Sj within NLO pQCD, who conclude that: i) there 
is no need to include intrinsic momentum effects, ii) the 
perturbation theory becomes unstable at lower kx when 
intrinsic momenta of partons is included, and iii) the data 
at lower energies are incompatible with those at higher 
energies, especially the p + Be data. Several other papers 
(see Ref. |54|) also discuss these aspects. 

On the other hand authors of Ref. [^5| have studied 
the effect of the so-called px broadening (Cronin effect) 
in proton-nucleus collisions. The LO pQCD is used along 
with a K factor and an intrinsic < kx > for the partons. 
The fact that the p + Be data is explained reasonably 
well by Vogelsang and Whalley Q using NLO pQCD, 
argued to have 'incompatible' normalization by Aurenche 
et al |34| l, and believed to require an intrinsic < kx > 
for partons as well as Cronin broadening of the intrinsic 
momenta of partons f35f| leaves the field open to diverse 
interpretations. From a purely empirical consideration 
also, it has been pointed out that |36| the lowest energy 
pp data for single photons are not consistent with the 
data at higher energies (and are too high). 

In the light of the above, we take the view that the 
estimates of the prompt photons given by Wong and 
Wang |l8fl , give the upper limit of these contributions. 

Summarizing, we find that the single photons mea- 
sured in the WA98 experiment seem to confirm the for- 
mation of quark gluon plasma in the collision and that 
most of thermal radiation at higher transverse momenta 
seems to come from the annihilation of quarks with scat- 
tering, which operates only if a plasma is formed. As 
expected, the slope of the spectrum provides a very good 
measure of the initial temperature reached in the colli- 
sion. 

This holds out the hope of a rich display of radiation 
of photons from the quark matter at RHIC and LHC en- 
ergies in collisions involving heavy nuclei, as much larger 
temperatures are likely to be attained there. The long 
life of the QGP phase at LHC energies will make it sen- 
sitive to such details like the transverse flow (within the 
QGP phase itself!), which will be of immense help in de- 
ciphering the properties of the quark-matter. 

Note added; further discussions 

We would like to take this opportunity to comment on 
some papers which have been posted 
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print archives after this paper was originally submitted. 
This discussion is necessary in view of the important con- 
clusions drawn in this work, which differs in detail with 
findings in these papers. 

The authors of Ref. |2{| have used a very simple model 
to parametrize the evolution of the plasma. A spherical 
(!) expansion of the plasma is envisaged which contin- 
ues to radiate photons during the entire life-time at a 
fixed (average, effective) temperature. While this may be 
useful to suggest that there is an additional production 



of photons, this approach is too simple to help us arrive 
at quantities like initial temperature, Tc etc. Moreover 
as, at the relevant nucleon-nucleon energy (^/s = 17.3 
GeV) for the WA98 data under consideration, there is 
no pp data, these authors further scale the predictions of 
PYTHIA for pp ( which required a K factor of 3.2 and 
intrinsic parton momentum for the E704 data at 19.4 
GeV) to the WA98 data for k T > 2.5 GeV. This fore- 
closes any hope to get information about the origin of 
these photons by assigning them to hard QCD interac- 
tions among partons whose distribution is obtained from 
structure function. 

The authors of Ref. ]3(| have used an early version of 
the transverse expansion code used in the present work, 
which was originally adopted from the work of Ruuska- 
nen and coworkers p3| . The model uses a energy-density 
profile which is unform upto the transverse radius R, a 
hadronic gas which consists of only 7r, p, w, r], and a\ 
mesons and nucleons, and uses the method of effective 
number of degrees of freedom at each temperature. This 
hadronic matter will have a much smaller number of de- 
grees of freedom at Tc, leading to a long lived mixed 
phase compared to the case of rich hadronic matter used 
in the present work. The overall life time of the sys- 
tem will then be larger, considerably enhancing the yields 
from the hadronic and the mixed phases. Thus one would 
need only a smaller contribution from the QGP phase to 
explain the data, as reported by these authors. It re- 
mains to be seen, how these results will behave when the 
a correction for the numerical factor of 4 jl2| in rates 
given by Ref. |ll| is made. Further, the method of tem- 
perature dependent number of degrees of freedom will 
lead to a speed of sound c s ~ l/v3 at all temperatures. 
A uniform energy density profile does not reflect the ac- 
tual situation either, as it would follow wounded-nucleon 
distribution used in the present work. 

Another aspect of this work is introduction of an initial 
transverse velocity. It is well known from the pioneering 
work of authors of Ref. |2^] that the < px > of the pro- 
duced particles can be arbitrarily increased if a strong 
initial transverse flow is assumed. A look at Fig. 3 of the 
present work also suggests that one can use a larger for- 
mation time (smaller initial temperature) and an initial 
flow to arbitrarily increase the large kx production. The 
authors of Ref. |3l]] have also introduced a large initial 
transverse velocity. 

However, it is known from the arguments of Ref. |25| l 
that in a head-on collisions of nuclei, introduction of ini- 
tial transverse velocity is not physically justified and one 
expects that, with the exception of the outer surface, 
the produced matter would be transversely at rest. At 
the most one may expect that the initial flow may be 
stronger near the surface, which expands against the vac- 
uum, there is no conceivable mechanism to provide a sig- 
nificant initial transverse collective motion, across the 
fluid. 

The initial scattering among partons will produce 
quarks and gluons pointing in random directions. In 



any given volume element their momenta would be uni- 
formly distributed in all directions, and rescatterings will 
then evolve a temperature and pressure. This tempera- 
ture and pressure (gradient) will initiate a flow when the 
plasma starts expanding against the vacuum. 

Both these works, Ref. [ |30|]3l] ], also include the LO 
pQCD predictions for the hard photons. While the au- 
thors of Ref. J|(| include the intrinsic &t as in the work 
of Ref. pH , they apparently they do not use a K- factor, 
though they use the LO pQCD. The authors of Ref. [§| 
use a iGfactor of 2 and find that they under-estimate the 
pp data at 19 GeV by a factor of 7 (implying an effective 
if- factor of ~ f4 (!) over the LO prediction). 

We have commented that the initial conditions de- 
duced here provide a quantitative description to the in- 
termediate mass dilepton spectra measured by the NA50 
group p0| . ft is of interest to understand the origin of 
the differences in the initial conditions inferred by us and 
those by the authors of Ref. p7\ , who report an initial 
temperature of about 200 MeV at an initial time of ~ 1 
fm/c. The fire-ball model used in Ref. jj7|, envisages a 
cylinder whose length and radius increase with time. The 
cylinder is assumed to be uniformly filled with plasma 
having temperature T{t). We know for sure that the 
profile of the energy density produced in such nuclear 
collisions can not be uniform, and this leads to additional 
gradients in the hydrodynamic evolution. The model also 
does not account for the fact that the speed of sound is 
large during the QGP phase, vanishingly small during 
the mixed-phase, and varying with temperature during 
the hadronic phase if a rich hadronic equation of state is 
used. Thus, for example, during the mixed phase their 
parameters a z and ax (which correspond to acceleration 
of the expanding surface) must vanish. Even though the 
parameters of the model are adjusted to give a transverse 
velocity equal to that deduced from particle spectra and 
a transverse size deduced from interferometry, it can not 
be expected to adequately reflect the rich history of evo- 
lution of the plasma formed in nuclear collisions, and by 
extension, the initial conditions. 

These simplifying assumptions provide that the con- 
tribution of the hadrons to radiations from the system is 
larger, necessitating only a small contribution from the 
QGP phase (and smaller initial temperature). 
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dN/dy=7^0 
t =1/3T =0.20 fm/c - 
T =335 MeV, T c =180 MeV: 
$WA98, Pb+Pb@SPS 1 
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FIG. 1. Single photon production in Pb + Pb collision 
at the CERN SPS. A chemically and thermally equilibrated 
quark-gluon plasma is assumed to be formed at to = l/3To 
which expands, cools, enters into a mixed phase and under- 
goes freeze-out from a hadronic phase. QM stands for radia- 
tions from the quark matter in the QGP phase and the mixed 
phase. HM, likewise denotes the radiation from the hadronic 
matter in the mixed phase and the hadronic phase. Prompt 
photons are estimated using NLO pQCD with the inclusion 
of intrinsic kr of partons (Wong and Wang [l8|). The (tail) 
ends of the arrows denote the upper limit of the production 
at 90% confidence limit. 
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T =0.20 fm/c, T =335 MeV- 
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FIG. 2. The sensitivity of single photon spectrum to crit- 
ical temperature. The solid curve is for Tc = 180 MeV, while 
the upper (lower) dashed curve is for 160 (200) MeV. 
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FIG. 3. The sensitivity of single photon (a) and pion 
spectrum J2(| (b) to initial time (temperature). The curves, 
from top to bottom, correspond to initial times of 0.2, 0.4, 
0.6, 0.8, and 1.0 fm/c for (a) and to 0.2 and 1.0 fm/c for (b). 



